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This  is  the  final  technical  report  covering  the  research  carried  out 
under  this  contract  from  March  1,  1975,  through  February  28,  1977.  The 
work  may  be  divided  into  these  four  categories: 

I.  Nonlinear  and  linear  stability  of  multistep  formulas, 

II.  Special  high-order  A-stable  methods, 

III.  Fractional  linear  difference  schemes  for  systems  of  quadratic 
differential  equations,  and 

IV.  Rotational  solutions  of  the  Josephson  phase  (Sine-Gordon-) 
equation. 

Here  we  summarize  our  results  obtained  in  these  areas.  A full 
account  is  given  in  the  research  papers  listed  in  Section  V and  referred 
to  hereafter  by  letters  in  brackets.  Some  of  these  papers  are  enclosed 
with  this  report  and  others  are  in  preparation.  Numbers  in  brackets  refer 
to  the  list  of  references  (Section  VI). 
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I.  NONLINEAR  AND  LINEAR  STABILITY  OF  MULTI  Sj’EP  FORMULAS 

/ 

A.  Nonlinear  Input-output  Stability  of  Conventional  Multistep  Formulas 
[a] 

If  a linear  multistep  formula  (LMF)  is  applied  to  a nonlinear  system 

i 

i 

i 

of  differential  equations,  y = f (y) , then  the  numerical  solution  { xn }, 
n = 0,1,...,  is  defined  by  a nonlinear  difference  equation,  * 

N(xn)  = 0;  (1) 

here  x approximates  y = y(t  ),  where  v(c)  is  an  exact  solution  of  the 
n n n 

differential  system,  t^  = nh,  li  > 0,  n = 0,1,...,  and 


N (x  ) = I a . x . 

n j=0  J n+J 


k 

-hr 

j=0 


e.f(x  ) 
1 n+j 


(2) 


If  by  d^  = -N(y^)  we  denote  the  local  truncation  error,  then  the  global 

error,  e = x -y  , satisfies  M(e  ) = d , where  M(e  ) = N(y  +e  )-N(v  ).  The 
nn  n n n n n n J n 

difference  equation  (1)  is  said  to  be  input-output  stable  if 

| |{en}|  | < K | | { d^ } | | in  some  appropriate  norm  ||*||,  where  K is  a constant 

independent  of  n.  In  order  for  the  concept  of  input-output  stability  to 
be  relevant  for  stiff  systems  of  differential  equations,  the  step-size  h 
must  be  thought  of  as  fixed  and  finite. 

We  have  proved  several  theorems  giving  sufficient  conditions  for 
input-output  stability  of  difference  equations  for  various  types  of  non- 
linearities,  notably  for  f's  satisfying  the  monoton i c i ty (dissipativi ty) 
condition 


< f (y+e) 


f(y). 


e > 


S P |e 


2 


(3) 


L, 


for  all  y and  e in  some  appropriate  set  and  for  some  u < 0;  here  <,>  is  an 


3 


i 


2 

arbitrary  scalar  product  and  |e)  = < e,e  >.  The  condition  (3)  implies 

stability  (dissipation)  of  the  differential  system.  We  have,  for  example, 
proved  the  following  result: 


i 

Theorem  1 : Consider  the  well-known  polynomials  p(t)  = I a , 

k . J-°  J 

a(0  = ^ & X3  associated  with  the  LMF  and  let  T denote  its  "root-locus 

j»0  3 

curve",  i.e.,  T = {q|  q = q(?)»  q(C)  = p(C)/a(0,  |c|  = 1).  Then,  if  the 


roots  Oj  of  a(?)  satisfy  |o^|  < 1,  j = l,...,k,  and  if  hp  < m,  where  m is 


the  deepest  incursion  of  T into  the  left  half  of  the  q-plane 

(m  = min  Req(O),  then  we  have  input-output  stability  both  in  the  l - and 

U 1=1 

in  the  maximum  norm. 

The  second  condition  of  Theorem  1 may  be  interpreted  by  saying  that 
the  differential  system  has  an  amount  of  dissipation,  measured  by  hp , 
sufficient  to  offset  the  lack  of  A-stability  [1]  of  the  LMF,  measured  by  m. 
Other  input-output  stability  theorems  have  been  obtained  for: 

• nonlinearities  f which  are  monotone  only  for  |e|  > B > 0, 

B sufficiently  large, 

• monotone  f's  with  sufficiently  large  negative  p, 

• nonlinearities  f which  are  gradients  of  concave,  scalar  functions 
(implying  that  the  Jacobian  matrix  3f/3y  is  symmetric  and  its 
spectrum  is  real), 

• nonlinearities  f satisfying  a Lipschitz  condition,  and 

• nonlinear  systems  composed  of  loosely  coupled  subsystems,  each 
of  which  falls  into  one  of  the  previously  mentioned  categories. 


A 


B.  Nonlinear  Stability  of  Formulas  with  Second  Derivatives  fb] 

The  nonlinear  input-output  stability  analysis  for  conventional  LMF 
described  in  the  preceding  paragraph  has  been  extended  to  the  LMF  involving 


second  derivatives: 


k 


2 a . x , , -h  2 3 . x , . + h 2 y . x , . = 0 . 

j=o  1 n+J  j=o  J n+i  j=o  -1  n+j 


(A) 


The  results  apply  to  stable  nonlinear  systems  x = f (x)  satisfying  the 

monotonicity  conditions 

h < f(x-t-e)  - f(x),  e > < | e | 2 , (p^  > 0) 

2 , . 2 
h < g (x+e)  - g(x),  e > 2 | o | , (p.,  > 0) 

< f (x+e)  - f (x) , g(x+e)  - g(x)  > < 0, 

where  g = f f and  f is  the  Jacobian  matrix  3f/3x. 
x x 

First,  simple  sufficient  conditions  were  derived  for  input-output 
stability  of  the  nonlinear  difference  operator  associated  with  the  dif- 
ferential system  and  the  integration  formula.  Rewrite  (A)  in  the  form 


2a. x - h 2 b . x ,+h  2c. k . = 0 

j=0  J n"J  j=o  J n_J  j=o  -1  n-J 


(5) 


where  a^  = «k_ j » b^  = 3^,  c^  = >k_j.  J “ 0,1,..., k,  n = k,  k+1,...,  and 

k i k 1 k \ 

let  r(z)  = 2 a . zJ  , s(z)  = 2 b.z  , and  t(z)  = 2 c.zJ  correspond  to  p(c), 

j=0  J j =0  J )=0  J 

, k i 

0(O , and  1(C)  = 2 Y,AJ»  respectively.  Then  the  above  mentioned  stability 

j = 0 ' 


condition  is  that 
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-1_  f uax  |£<£il . A|\  + i f „ax  |LSii 
a+bV"2  \ 1*1-1  Wl>  J B Vui-1  W1) 


B I ) < 1,  (6) 


where  A and  B are  free  positive  parameters  which  are  at  one's  disposal 
and  which  are  to  be  suitably  chosen. 

Second,  the  criterion  (6)  was  applied  to  the  A-stable  formulas  of 
Enright  [2]  with  k=l,  p=3  and  with  k=2,  p=4 , respectively.  A computer 
search  procedure  was  used  to  find  values  of  A,B  for  which  (6)  is  satisfied 


for  various  amounts  of  dissipation,  i.e.,  various  values  of  and  p9.  For 
example,  Enright's  formula  with  k=l,  p=3, 

* - , , - jj  (2i  + i - 0,  (7) 


was  found  to  be  input-output  stable  with  a moderate  amount  of  dissipation; 

M1  " W2  = 1*3* 


C . Liapunov-Stability  of  Conventional  Multistep  Formulas  [ c ] 


If  {x  } and  {x  +z  } are  any  two  solutions  of  the  difference  equation 
n n n 

(1),  then  {z^}  satisfies  the  variational  difference  equation 

N(x  ) - N(x  +z  ) = 0.  The  LMF  giving  rise  to  this  variational  equation  is 
n n n 

said  to  be  Liapunov-stable  if,  in  some  appropriate  norm  ||*||  defined  by  a 

Liapunov  function,  | |{z  I 1}  < “.  It  is  said  to  be  (G.p)-stable  in  the  sense 

n 

of  Dahlquist  if  for  any  f satisfying  the  monotonicity  condition  (3)  it  is 
Liapunov-stable  with  respect  to  a norm  defined  by  a positive  definite  quad- 


ratic vector  form  Z g,,  < z , z , > , where  G«(g,,)  is  a positi 

j )=0  ij  n+k-i*  n+k-j  Bij 
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definite  real  symmetric  kxk-matrix.  The  LMF  is  said  to  be  G-stable  if  it  is 
(G, 0)-stable . A-stability  is  necessary  for  its  non-linear  analogue, 
G-stability . 

By  a simple,  explicit  construction  of  a quadratic  Liapunov  function 
we  proved  the  following  results: 

Theorem  2:  Any  member  of  the  four-parameter  family  of  all  three-step 

(k  = 3)  LMF  which  are  second-order  accurate  (p  = 2)  is  G-stable  if  and  only 
if  it  is  A-stable. 

• Any  three-step  (k  = 3),  third-order  (p  = 3)  LMF  whose  polynomial 
a(?)  satisfies  assumption  1 of  Theorem  1,  as  well  as  certain 
other,  easily  verifiable  algebraic  constraints  between  its  co- 
efficients, is  (G,u)-stable  for  any  p below  an  explicitly  given, 
negative  bound  which  depends  on  the  choice  of  the  formula.  For 
example,  the  formula 

-lOx  + 42x  -78x  _+46x  ^-h(f  -f  ,.-f  _u0+25f  = 0, 

n n+1  n+2  n+3  n n+1  n+2  n+3 

a "neighbor"  of  the  three-step  backward  differentiation  formula, 

is  (G.p)-stable  for  p <-36/31. 

A survey  of  work  on  the  stability  (mostly  linear)  and  accuracy  of 
numerical  methods  for  stiff  differential  equation  (d]  was  presented  by 
the  principal  investigator  at  the  workshop  on  stiff  differential  equations 
held  at  the  Air  Force  Weapons  Laboratory  at  Kirtland  Air  Force  Base,  N.M., 
on  May  6 and  7,  1976. 

A survey  of  recent  work  on  nonlinear  stability  of  integration  methods 
is  in  preparation  [ e ] . 
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II.  SPECIAL  HIGH-ORDER  A-STABLE  METHODS 

A.  Averaging  [f] 

High-order,  A-stable  numerical  solutions  to  stiff  systems  of  differ- 
ential equations  may  be  generated  by  a)  calculating  several  A-stable  solu- 
tions of  order  of  accuracy  one  or  two,  each  of  which  is  produced  by  a parti- 
cular member  of  a family  of  LMF  depending  on  instrinsic  parameters,  and  b) 
forming  a suitable  "weighted  average"  (linear  combination)  of  these  solu- 
tions. The  recipe  for  forming  high-order  averages  of  low-order  LMF  solu- 
tions is  found  by  considering  an  asymptotic  expansion  of  the  global  trun- 
cation error  as  a function  of  the  formula  parameters. 

• An  efficient  procedure  was  developed  for  the  systematic  compu- 
tation, to  high  orders  in  h,  of  an  asymptotic  expansion  of  the 
global  truncation  error  of  an  LMF.  This  technique  was  applied 
for  finding,  to  O(h^),  the  error  expansion  for  the  four-parameter 
class  of  all  three-step,  second-order  LMF.  From  a two-parameter 
subclass  of  this  class,  A-stable  (in  fact,  G-stable)  methods 
of  orders  4,5,  and  6 have  been  derived  which  require  the  aver- 
ging  of  as  few  as  2,  3,  and  4 LMF  solutions,  respectively.  The«e 
methods  have  been  programmed  in  APL  and  successfully  tested  on 
linear  and  mildly  non-linear  stiff  systems.  The  success  of  the 
numerical  tests  hinged  partly  on  the  use  of  a suitable  starting 
procedure  we  developed  for  the  step-by-step  solution  of  the  dif- 
ference equations. 

B . A-stable  Integration  Formulas  with  Second  Derivatives  [ g ] 


Whereas  for  conventional  linear  multistep  formulas  (LMF)  A-stability 
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I 


is  incompatible  with  orders  of  accuracy  p > 2 [1],  there  do  exist  A~stab]e 
LMF  involving  second  derivatives,  i.e.,  formulas  of  the  type  (4)  of  as 
high  an  order  of  accuracy  as  four  [2,3].  A criterion  was  developed  for 
testing  A-stability  of  formulas  of  class  (1),  which  also  served  as  a basis 
for  the  ci  priori  construction  of  a special  class  of  A-stable  formulas  of 

k 

this  type.  This  criterion  states  that  if  i)  the  formula  (4)  is  A^-stable  , 

k 

i.e.,  the  roots  x.  of  the  polynomial  x(c)  = £ satisfy 

1 j=0  3 

| t ^ | < 1,  i=l,...,k;  and  ii)  under  the  two-valued  map  r,  -*■  q,  defined  by  the 
characteristic  equation 

p(4)  - qo(4)  + q^xU)  = 0,  (8) 

the  image  set  T of  the  unit-circle  of  the  £-plane  (sometimes  referred  to 

as  the  root-locus  curve)  satisfies  r c {q|Req  > 0},  then  the  formula  (4) 

is  A-stable.  This  criterion  is  based  on  an  earlier  result  of  ours  [4]  and 

generalizes  a criterion  valid  for  conventional  LMF  [5],  i.e.,  formulas  with 

x(0  i 0.  It  can  be  implemented  efficiently  by  using  the  algorithms  of 

Routh  [6]  and  those  described  by  Duffin  [7]. 

The  criterion  mentioned  above  has  been  used  to  a priori  construct 

the  three-parameter  family  of  all  A-stable  two-step  (k  = 2)  formulas  of 

type  (1)^  which  have  p = 4.  The  formulas  of  Enright  [2]  and  those  proposed 
by  Jackson  and  Kenue  [8]  are  members  of  this  family. 


The  present  definition  of  A^  -stability  is  analogous  to  that  for  con- 
ventional LMF  given  in  [4]. 


t 


Except  possibly  for  some  marginally  A-stable  ones. 


[g] 

If  the  formula  (1)  is  applied  to  the  test  equation 


x = Ax,  A = const.  (9) 

the  fundamental  solutions  of  the  resulting  difference  equation  are  of  the 
form  {?n,  n = 0,1,2,...}  where  C is  any  one  of  the  solutions  of  the  character- 
I istic  equation  (8)  with  q = hA , provided  these  roots  are  distinct.  If  one 

of  them  satisfies  C = e , i.e.,  if 

p(0  - log;  o(0  + (log;)2rU)  = (10) 

then,  for  this  t,  f^11}  is  an  exact  discrete  solution  of  (9),  i.e., 

Cn  = e^U  = x(t  ),  t = nh,  where  x(t)  = e^t  is  the  exact  fundamental  solu- 
n n 

tion  of  the  differential  equation  (9).  Then,  as  far  as  this  "principal" 
root  is  concerned,  the  formula  (4)  is  exponentially  fitted  [9,  10]  at 
q = log  C. 

Let  R(0  denote  the  left  side  of  (10).  For  any  given  |R(c)|, 
is  a measure  of  how  far  the  formula  is  from  being  exponentially  fitted. 

1 2 

Thus,  M = / [R(01  ds  is  an  i?-measure  of  global  accuracy  of  he  formula 

o 

(4)  with  respect  to  the  family  of  test  problems  x = Ax,  S A S 0.  The 
measure  M can  be  computed  explicitly  as  a function  of  the  formula  parameters. 

For  the  particular  three-parameter  class  of  formulas  defined  under 
II. B the  following  problem  was  solved:  Optimize  the  formula  with  respect 

to  M over  all  three  parameters,  subject  to  the  A-stability  constraints.  The 
unique  optimal  formula  was  found  to  be 


i 
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1236x  -241620x  J,+240384x  ^,,-h(-618x  +119574x  , ,+120192x  _) 
n n+1  n+2  n n+1  n+2 

+h2(103x  -20135x  L1+20032x  = 0. 

n n+1  n+2 


III.  FRACTIONAL  LINEAR  DIFFERENCE  SCHEMES  FOR  SYSTEMS  OF  QUADRATIC 
DIFFERENTIAL  EQUATIONS  [h,i,j] 

We  are  concerned  with  the  numerical  solution  of  quadratic  systems 

of  ordinary  differential  equations 

.i  i a 8 i a i 
x = a „x  x + b x +c, 

ocg  a 

by  quadratic  systems  of  difference  schemes 

, i 6/x,  i.a.,,.  i a , . g . . , i a,  . , i 

(-r  _u  (t)  + s )u  (t+h)  = t „u  (t)u  (t)  + u u (t)  + v . 
ag  a ag  a 

Here,  a and  6 are  summed  from  1 to  n,  and  i ranges  from  1 to  n. 

In  vector  notation  the  initial  value  problems  for  these  equations  art 


x = x*x  + Bx  + c,  x(0)  = £, 


(ID 


(-R[u(t ) ] + S) u ( t+h)  = T[u(t)]u(t)  + Uu ( t ) + v , u (0)  = £.(12) 

In  the  above,  x,  c,  £,  u and  v are  vectors  in  Rn;  B,  S and  U are  n«n 

matrices;  * is  multiplication  in  the  algebra  <!?(*)  = (a!,);  and  R[u]  and 

J k 

T[u]  are  n»n  matrices  which  are  linear  functions  of  u. 

Wheh  T[u]  = 0,  the  scheme  (12)  can  be  written  as 


u(t+h)  = (-R [ u ( t ) ] + S)  J(Uu(t)  + v), 


(13) 


which  we  call  a fractional  1 inear  scheme.  We  have  obtained  the  following 
results : 

1)  If  the  scheme  (13)  converges  as  h -*■  0 for  all  £,  then  it  must 
converge  to  a quadratic  system  (11). 


11 


2)  For  every  quadratic  system  (11),  there  exists  a fractional  linear 
scheme  (13)  which  approximates  (11)  to  second  order  and  which  converges  to 
(11)  as  h -*■  0. 

3)  If  the  coefficients  of  the  scheme  (13)  satisfy 


R(h) [u ] p = / V (x) ((V(t)u)*(V(T)p))dr , (14) 

o ' i 

S (h)  = V(h)  (15) 

U (h)  = -R(h) [V (h) z (h) ] + I,  (lb) 

v (h)  = V (h) z (h) , (17) 

where  z = z(t)  and  V = V(t)  are  defined  by 
z = z*z  + Bz  + c,  z(0)  = 0, 


V = -V (B+2R* [z ] ) , V(0)  = I,  (R*[z]x  = z*x)  , 

then  (13)  is  a second  order  accurate  convergent  scheme  for  (11). 

4)  The  scheme  (13)  is  the  exact  scheme  for  (11)  if  and  only  it  the 
scheme  coefficients  satisfy  (14)-(17)  and  the  system  (11)  is  a P-system. 
That  is,  the  solution  of  (11)  is  given  by 

x(t,0  = (-R(t)mv(t)+V(t))_1((-R(t)[V(t)z(t)]+m  + V(t)z(t)). 

5)  If  the  coefficients  of  (13)  satisfy  (14)—  (17)  to  second  order 

in  h,  then  (13)  is  also  a second  order  accurate  convergent  scheme  for  (11). 

6)  The  matrix  Riccati  equation 

X = XAX  + BX  + XC  + D, 

which  is  of  great  practical  importance.  Is  a P-system  and  therefore  has 
a fractional  linear  exact  scheme.  Further,  this  exact  scheme  can  be  con- 


structed by  solving  a 1 lnear  system  of  ordinary  differential  equations 
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on  the  Interval  [0,h]. 

These  results  give  a rather  complete  analytic  theory  of  fractional 
linear  difference  schemes  for  quadratic  systems.  To  test  the  computational 
effectiveness  of  these  schemes,  two  computer  codes  are  being  developed  im- 
plementing these  schemes:  one  for  two-dimensional  systems,  and  one  for 

arbitrary  matrix  Riccati  equations. 


IV.  ROTATIONAL  SOLUTIONS  OF  THE  JOSEPHSON  PHASE  (STNE-GORDON)  EQUATION [k  ,2  ] 
A.  Problem  Formulation 

A Josephson  junction  consists  of  two  superconductors  separated  by 
an  extremely  thin  dielectric  barrier.  The  order  parameter  <j>  (phase  dif- 
ference between  the  wave  functions  of  the  two  superconductors)  satisfies 
the  Josephson  equations  [11].  In  a mathematical  limit  situation  where  <j> 
depends  only  on  one  space  dimension  and  on  time,  Josephson's  fundamental 
equations  combined  with  Maxwell's  equations  give  rise  to  a nonlinear  damped 
wave  equation  for  <f>  referred  to  as  the  Josephson  phase  equation  or  the 
"sine-Gordon"  equation.  In  dimensionless  form  this  equation  is 

A«f>  i | ^ -a  — = k sin  ij> . (18) 

3x“  3t“ 

For  the  purpose  of  describing  the  Josephson  ac-effect,  one  is  seeking 


rotational  solutions  of  (18);  i.e.,  solutions  of  the  form  <f  = wt  + ip , where 
ip  is  periodic  in  t of  period  (2ir/w),  for  some  w . 


The  dimensionless  variables  are  those  introduced  in  [12],  except  that 
distances  are  measured  in  units  of  the  I./n  where  1.  is  the  junction  length, 
rather  than  in  units  of  the  Josephson  penetration  depth  \ . 
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In  a situation  where  a constant  voltage  is  applied  to  one  end  of 
the  junction  (x=0) , the  equation  (18)  is  to  be  solved  subject  to  the 
boundary  conditions  [12] 


x=0 


wt 


u 


V 

e 


Mi 


X=TT 


where  denotes  an  externally  applied  magnetic  field.  This  case  is  re- 

2 2 

ferred  to  as  the  "voltage-driven  case."  Note  that  k = L /(ttX  ) where  X 

J J 

is  the  Josephson  penetration  depth,  and  the  o of  [12]  equals  (nA^/L)  times 
the  present  a.  The  total  current  I drawn  by  the  Junction  is  the  time  average 
(over  a period)  of  the  ac  total  current  I(t)  defined  by  the  additional  con- 
straint 


| - F I n = KO. 

3x  'x=ir  3x  'x=0 


3 <J> 


The  nonlinear  current-voltage  characteristic  of  the  junction  is  then  given 
by  I = I (u>)  . 


If  the  junction  is  driven  by  a constant  total  imposed  current  Ie> 
one  is  looking  for  rotational  solutions  of  (18)  with  unknown  w satisfying 
the  boundary  conditions 


M I 

3x  ' x-0 


-(He  + V’ 


^ | = -H  • 

3x  1 x-JT  e 


Then  the  voltage  V across  the  junction  at  x » 0 is  the  time  average  of  the 
ac  voltage  V(t)  defined  by  V(t)  * (3<^>/3t)  Ix_q-  This  case  is  referred  to  as 

the  "current-driven  case." 


r 


I 

s 

J 

■ 

* 

f 


B.  Perturbation  Solutions 

a)  Analytic  solutions  valid  for  small  nonlinearities  (i.e.,  small 
values  of  k,  corresponding  to  junction  lengths  L which  are  small  compared 
to  Aj)  were  obtained  by  perturbation  methods  in  both  the  voltage  and  the 
current  driven  cases.  They  exhibit  the  typical  "resonance"  behavior  of  the 
I-V-characteristics  observed  experimentally  and  predicted  numerically  [14]. 

b)  A constructive  perturbation  procedure,  similar  to  those  used  in 
bifurcation  theory  and  valid  in  the  limit  of  strong  dissipation  (o  -*•“>), 
was  developed  for  the  voltage  driven  case. 

C.  Existence,  Uniqueness,  and  Stability  Results 

If  one  writes  <t>  = <Pq  + ip,  where 

2 

= wt  - kx  + px“ 

satisfies  the  formal  limit  equation  A4> ^ = 0 as  tc  -*•  0,  and  where  k = ou)  + H , 
p = ou)/2,  and  ip  is  periodic  in  t,  then  in  the  voltage  driven  case  ip  satis- 
fies 


Aip  = k sin(i[>Q  + ip)  , 


I = — I 

1 x=0  3x  1 x=ir 


0. 


(19) 


The  following  results  were  proved  in  this  case: 

a)  For  any  a?4  0 there  exist  (2ir/u))-periodic  solutions  of  the  problem 
(19).  Hence,  for  o^O,  there  exist  (possibly  multiply  branched  [13]) 
I-V-characteristics  for  the  Josephson  junction. 

b)  For  moderate  amounts  of  dissipation  relative  to  the  strength  ol 


the  nonlinearity  (i.e.,  moderately  large  values  of  o relative  to  k),  the 
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periodic  solution  of  (19)  is  unique  and  globally  asymptotically  stable. 
Hence  solutions  with  arbitrary  initial  data  tend  - exponentially  in  t - 
to  <)>  = (j>n  + t|)  with  periodic.  For  example,  if  < = 5 (corresponding  to 

U O 

L/Xj~. 9)  then  a t 9/16  is  sufficient  for  uniqueness;  similarly,  with 
2 1 

k = 1 / 4ir  (corresponding  to  L/X^  = uniqueness  is  assured  for  0 £ .06. 
These  results  validate  the  perturbation  procedure  discussed  under  IV. h. a 


in  the  limit  k -*  0. 
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